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We study the effects of a power-law trapping potential on the scaling behaviour of 
the entanglement at the quantum critical point of one-dimensional (ID) lattice particle 
systems. We compute bipartite von Neumann and Renyi entropies in the presence of 

"j^ , the trap, and analyze their scaling behaviour with increasing the trap size. 

'tr. ' As a theoretical laboratory, we consider the quantum XY chain in an external 

transverse field acting as a trap for the spinless fermions of its quadratic Hamiltonian 
representation. We then investigate confined particle systems described by the ID 
Bose-Hubbard model in the superfluid phase (around the center of the trap). 

i-rt ■ In both cases conformal field theory predicts logarithmically divergent bipartite 

^ ' entanglement entropies for the homogeneous systems without trap. The presence of 

the trapping potential breaks conformal invariance, affecting the critical behaviour of 
the homogeneous system. Our results show that the bipartite entanglement entropies 

^^ diverge logarithmically with increasing the trap size, and present notable scaling 

^ . behaviours in the trap-size scaling limit. 

o 
in 

. • ' 1. Introduction 

in 

o, 

O . The spatial entanglement of one-dimensional (ID) lattice systems near their quantum 

critical point is a physical issue which has attracted much interest, see, e.g., Ref. [1]. 
Many studies have been devoted to quantify it, by computing von Neumann or Renyi 
entanglement entropies of the reduced density matrix of a subsystem. In a homogeneous 
system of length L whose quantum critical behaviour is described by a two-dimensional 
(2D) conformal field theory (CFT), the bipartite entanglement entropy S{Ia, L) increases 
logarithmically [21 |3l IH |5] with increasing the size Ia of the subsystem at the critical 
point, i.e., S ~ In/^, or with increasing the spatial length scale ^ when approaching 
the critical point, i.e., S" ~ In^ (if 1 <^ ^ ^ Ia)- In both cases the coefficients of the 
logarithms are related to the central charge of the corresponding CFT. 

Quantum lattice particle systems are realized in experiments of cold atoms in 
optical lattices of any effective spatial dimension d < 3, which allow to investigate the 
interesting interplay between quantum and statistical behaviours in particle systems, 
see, e.g., Refs. P[7llHll9l[ini[IIl[12l[13l[Tlll5l[l6l[I7]. An important feature of these 
experiments is the presence of a confining potential which traps the particles within 
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a limited spatial region of the optical lattice created by laser-induced standing waves. 
The theoretical framework |T8] is based on the Bose-Hubbard (BH) model [T9] in the 
presence of a confining potential coupled to the particle density, i.e., 

i^BH = --5^(&l&j+&J&i)+- J]ni(ni-l)+/i^ni + ^V(r,K, (1) 

{ij) i i i 

where (ij) is the set of nearest-neighbour sites, Ui = h\hi is the particle density operator. 
We consider a power-law trapping potential 

V{r) = v^rP = {r/iy, (2) 

where r = \x\, v and p are positive constants and / = 1/v is the trap size. Experiments 
are usually set up with a harmonic potential, i.e., p = 2. Far from the origin the potential 
V{r) diverges, therefore (nj) vanishes and the particles are trapped. The inhomogeneity 
due to the trapping potential strongly affects the phenomenology of quantum transitions 
in homogeneous systems. The correlation functions develop critical behaviours only in 
the limit of large trap size. This critical behaviour can be described in the framework 
of the trap-size scaling (TSS) theory pUf fIT\. 

In this paper we focus on the quantum entanglement in lattice particle systems 
confined by an external potential, such as the BH model ([1]). We investigate how the 
scaling behaviour of the entanglement entropies determined by CFT changes in the 
presence of the inhomogeneity induced by the confining potential. 

As a theoretical laboratory we consider the quantum XY chain in an external space- 
dependent transverse field, acting as a trap for the spinless fermions of its quadratic 
Hamiltonian representation. This model presents a quantum critical behaviour 
belonging to the 2D Ising universality class, thus described by a CFT with central 
charge c = 1/2. In the presence of the trap and for large trap sizes, the quantum critical 
behaviour shows TSS [Sni El]. The presence of the trap induces a length scale ^ at the 
critical point, which scales nontrivially with increasing the trap size, i.e., as ^ ~ /^ with 
a trap exponent 9 = p/{p + 1). 

We then consider the ID BH model ([1]) in its hard-core limit, i.e., ?7 — )• oo, and study 
the bipartite entanglement entropy in the gapless superfluid phase, whose continuum 
limit is described by a 2D massless bosonic field theory, thus a CFT with c = 1. In 
this case the trap-size dependence is made more comphcated by new features [22]: (i) 
a multiscaling phenomenon, i.e., the existence of different length scales which diverge 
with distinct power laws in the TSS; (ii) the presence of level crossings of the lowest 
states at finite trap size, which gives rise to a peculiar modulation of the amplitudes of 
the asymptotic power-law behaviours. 

We analyze the scaling behaviour of bipartite von Neumann and Renyi 
entanglement entropies in these ID lattice models. We consider chains with even sites 
L and open boundary conditions, and a trap of size / centered between the middle sites 
of the chain. We divide the chain in two connected parts of length I a and L — Ia, and 
consider their von Neumann entropy 

Si{Ia]L) = Si{L -Ia]L) = -TilpAlnpA] (3) 
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and their Renyi entropies 

S^{Ia]L) = S^{L-Ia;L) = -^InTrp^ (4) 

1 — a 

where pA is the reduced density matrix of one of the two subsystems. Renyi entropies are 

useful because they provide information on the spectrum of pA [23] . The von Neumann 

entropy ([3]) is obtained by taking the hmit a — )■ 1 in (jlj). In the presence of a trap of 

size I, the bipartite entanglement entropies of subsystems of size I a = L/2 — x have a 

finite L — )■ oo limit, depending on x and /. We study their scaling behaviour in the 

large-/ limit. Note that the p ^^ oo limit of the confining potential becomes equivalent 

to a homogeneous chain of size L = 21 with open boundary conditions (more precisely, 

to a chain with even L = 2[l\ sites). Thus, in the p — )■ oo limit we must recover the 

entanglement entropies of homogeneous systems with open boundary conditions, which 

are determined by CFT. 

In our numerical study we exploit the quadratic spinless fermion representations of 
the XY chain and the ID hard-core BH model, which allow us to perform computations 
for very large systems, since they only require the diagonalization of a L x L matrix 
where L is the number of lattice sites. We compute the entanglement entropies using the 
method outlined in Refs. [2U |25] , which exploits the quadratic fermionic representation 
and can be also applied in the presence of the trapping potential. We obtain numerical 
results for chains of size L and open boundary conditions, with a trap of size / centered 
between the middle sites of the chain; we choose L large enough to have negligible 
finite-L effects; we are able to obtain very accurate results for / up to O(IO^). 

Our results show that the bipartite entanglement entropies at the critical point 
diverge logarithmically with increasing the trap size. Moreover, they present notable 
scaling behaviours in the trap-size scaling limit. 

We mention that the bipartite entanglement entropies in the XY and XX model with 
gradients, i.e., in the presence of a linear external field, have been recently investigated 
in Ref. [26]. 

The paper is organized as follows. In Sec. [2] we presents results for the quantum 
XY chain in an external space-dependent transverse field. Sec. [3] is dedicated to the ID 
hard-core BH model. In Sec. Hlwe draw our conclusions. 

2. Bipartite entanglement in the XY chain with a space-dependent 
transverse field 

The quantum XY chain in a transverse field is a standard theoretical laboratory for 
issues related to quantum transitions, see, e.g., Ref. [27]. An inhomogeneity analogous 
to the one arising from a trapping potential in particle systems can be achieved by 
considering a space-dependent transverse field, i.e., 

H^y = -H ^[(1 + ^^^^^t,i + (1 - 7)^^+1] - f^crt - VixMt, (5) 
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where a* are the Pauh matrices, 7 7^ and V{x) = (|x|//)^, cf. ([2]). The XY chain 
can be mapped into a quadratic Hamiltonian of spinless fermions by a Jordan- Wigner 
transformation, obtaining 

H = Y.^clA^''^ + \{c\B,,c] + h.c.)], 

Aij = 26ij - 5i+ij - 6ij+i + 2Q{xi)6ij, /g\ 

Bij = —7 ((^i+ij — Si,j+i) , 

Q{x) = Ji + V{x), fL = fi-l. 

In this picture fi plays the role of chemical potential for the fermion c-particles, and the 
space-dependent field V{x) acts as a trap for the c-particles, making their local density 
(ui) = {clci) vanish at large distance. In the following, / will be considered as the trap 
size. 

In the absence of the trap, the model undergoes a quantum transition at /i = 
yU — 1 = in the 2D Ising universality class, separating a quantum paramagnetic phase 
for p, > from a quantum ferromagnetic phase for /2 < 0. Around p. = 0, the quantum 
critical behaviour shows a diverging length scale, ^ ~ l/^l '^j ^iid a vanishing energy 
scale, A ~ ^~^, where z and z/ are universal critical exponents: z = 1 and u = 1/y^ = 1 
{y^ is the RG dimension of /x) . 

The presence of the external field V{x) gives rise to a space inhomogeneity, thus 
affecting the critical behaviour of the homogeneous system, which is only recovered in 
the limit of large trap size. In this limit the critical behaviour can be described in the 
framework of the trap-size scaling (TSS) theory [20]. The TSS Ansatz for the asymptotic 
behaviour of the free-energy density in the presence of a confining potential ([2]) is 

F(/i, T, /, x) = l-'^^+'^T{jjl'/\ Tf', xl-'), (7) 

where x is the distance from the middle of the trap, and 6 is the trap-size exponent which 
determines how the length scale of the critical modes at the critical point diverges with 
increasing trap size, i.e., ^ ~ l^. RG arguments and analytic calculations show that [20] 

e = p/{p+i). (8) 

The TSS limit can be analytically derived within the quadratic spinless fermion 
representation ([6]), obtaining a Schrodinger-like equation for the lowest states, after 
introducing the rescaled spatial variable X = 7~^/^Z~^x, and the rescaled deviation 
from the critical point fir = 7~^/^/i (see Ref. [20] for details). The 7 dependence gets 
absorbed by such rescalings leading to a universal TSS limit. Therefore at the quantum 
critical point /2 = and T = 0, any length scale related to the critical modes must have 
the asymptotic behaviour 

^ = a7^V [1 + 0(/-')] (9) 

where the amplitude a depends on the power p of the potential, but not on 7. The 
power law of the scaling corrections is inferred from the analysis of the corrections to 
the continuum Schrodinger-like equation describing the TSS limit. 
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We want to study the quantum entanglement of the ground-state in the presence 
of the external space-dependent potential (|2]). For this purpose we analyze the trap-size 
dependence of the bipartite von Neumann and Renyi entanglement entropies, defined 
in equations ([3]) and (jl]), in the framework of the TSS theory. 

Before presenting the results in the presence of the trap, let us recall some 
established results for the homogeneous XY chain with open boundary conditions, which 
represents the formal p — )■ oo limit of the model with the confining potential. They can 
be obtained by exploiting the CFT that describes its continuum limit. Setting 

Ia = L/2-x (10) 

(we consider a chain with an even number of sites L), the Renyi entropies at the critical 
point /i = behave as [1] 



In L + In cos ( — ) + e^ 
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where 



a = ci±|^ (12) 

and c is the central charge: c = 1/2 for the 2D Ising universality class. The constant e^ 
can be derived using the results of Refs. [U |2S1 [2H1 ED] , obtaining 

e„ = ln7 + ln(8/7r) + y„, 

y^= r^ \-^- (^ U] -1 e-H . ^''^ 

Jo ^ 1 — a^ \asinht/a sinhty sinht 
The behaviour of the von Neumann entropy can be obtained by taking the limit a — )■ 1 
of ( ITTl) . Corrections to the asymptotic behaviour (ITTl) are expected to be 0(L^^/") 
at fixed x/L [30l [31] . They become logarithmic for a — )■ oo. In the case of the von 
Neumann entropy, i.e., a — )■ 1, corrections are 0(L^^)|j 

In the presence of the trap of size /, the dependence on L rapidly disappears for 
L sufficiently larger than the trap size /. We define entanglement length scales C,a from 
the half-lattice entanglement entropies. We write (at p. = 0) 

^„,i/2 = lim S^{L/2; L) = C, (In^a + e„) , (14) 

where Cq, is the constant defined in (TT3l) . We identify the argument of the logarithm as 
the entanglement length scale ^a, which will generally depend on a and on the power p 
of the potential (I2l). However, by construction, their p — )■ oo limit is ^^ — ^ 2/ = L for 
any a. Since the entanglement length scales are clearly related to the critical modes, 
they are expected to have the asymptotic behaviour ([9]) at the critical point, although 
their scaling corrections may turn out to be different, depending on a. 

The approach to their asymptotic behaviour ^q, ~ l^ may be guessed by recalling 
that the p ^ oo limit is equivalent to a homogeneous chain of size L = 21 with open 
boundary conditions, where the corrections to the asymptotic behaviour fITT]) of the 



I It is worth mentioning that, using the relation [29] Sa{lA',L)\j^i = ^Sa{'2lA','2L)\^=Q and results 
from Refs. [21 [H], we obtain S'i(L/2; L) = ^ [inL + d + (2 - 37r)/(4L) + 0(l/i2)] for 7 = 1. 
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Figure 1. The ratio £,i/i"f^^^l^) vs l^^ for p = 2, 4, several values of 7, and trap sizes 
in the range 10^ < / < 10^. The lines are drawn to guide the eyes. The crosses along 
the y-axis represent the 7-independcnt I — )• 00 extrapolations, see text. 



Renyi entanglement entropies are [31] 0(L~^/"). Then, the correspondence L ^ ^ r^ l^ 
suggests the presence of 0(^~^/'^), thus 0{l~^^°'), corrections. In the case of the length 
scale derived from the von Neumann entropy, i.e., for a — )■ 1, this argument leads to the 
same scaling corrections of equation i^, which are expected for any length scale defined 
from standard correlation functions, such as Gij = {a^a^). 

The scaling (jQ]) is fully supported by numerical results for the half-lattice von 
Neumann entropy. |§| Figure [1] shows the results of ^1 for p = 2 and p = A, several 
values of 7, and for trap sizes in the range 10^ < / < 10^. In order to estimate the 
amplitude ai of the asymptotic behaviour 

^, = a,Y^^l'[l + h^l-' + ...], (15) 

we fit the data of the ratio ^1/(7^^^/^) to the Ansatz ai + 61,1/^^ + 61, 2^^^^, which turns 
out to be optimal, and check the stability of the results by increasing the minimum 
value of / allowed in the fit. The analysis of the p = 2 data gives ai ^ 1.97431 
for all 7 = 6,4, 2, 1, 1/2IIJ confirming its independence of 7, and 611 = 61,17^/^ with 
611 ~ —2.677 independent of 7. For p = 4 we obtain ai ^ 2.18316 for both 7 = 4, 1. 

In the case of the entanglement length scales C,a defined from the Renyi entropies, 
we generally expect ^q, ~ l^ asymptotically, with 0(/~^/") corrections. Assuming 
universality of the TSS limit, we expect that the ratios C,a/^i approach universal values 



a/ei = Ra + 0(/ 



-o/a\ 



(16) 



§ The effects of the finite size L of the chain are under complete control in our numerical calculations 
at fixed trap size. This is easily achieved by comparing data at fixed trap size for increasing chain 
size, due to the fast convergence to the infinite-chain limit, see also later. We choose L large enough 
to have negligible finite-L effects, so that the results that we present can be effectively considered as 
infinite-chain results with great accuracy. Therefore, the systematic error related to the finite size of 
the chain is totally negligible in our analyses. 

II Here and in the following the uncertainty on the numerical estimates is at most one on the last 
reported figure. 
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Figure 2. Ratios ^q/^i of entanglement length scales vs. Z ^/"a in the presence of an 
external potential with p = 2, 4, for a = 2, 3, several values of 7, and trap sizes in the 
range 10^ < Z < 10''. 



and therefore 

ia = a^^"n' [I + 6„,i(7)r^/" + 6„,2(7)^"''/" + •••] , (17) 

with Oq dependent on p but not on 7. Note that i?^ — )■ 1 and Oq — )■ 2 for p — )■ 00. 

This scenario is confirmed by numerical resuhs for the ratios E,a/^ii see figure [2] 
where we show resuhs for p = 2, 4 and a = 2, 3. They fit the general behaviour 

iJii = R^ + Y,c,n"^, (18) 

i=i 
and provide consistent values of R^ for different values of 7, supporting universality. 
For example, we find R2 ^ 0.9889, R^ ^ 0.987 for p = 2, and i?2 ~ 0.9970 and 
i?3 ~ 0.996 for p = 4, which are slightly smaller than one, and appear to approach one 
with increasing pl^fl 

Moreover, the results for p = 2 indicate that the amplitudes 6a,i of the leading 
0(/^^/") scaling corrections in equation f lT7|) behave as ha^i = &a,i7~''''" where k = 4/3 
and ba,i is independent of 7, with great numerical accuracy. Analogous results are found 
for p = 4, but with k = 6/5. Taking into account that ba,i ~ 7^^/°" for p —^ 00 (as 
inferred by numerical computations), we may guess k, = {p + 2)/{p + 1). Therefore, we 
have that the ratios 6^ i/&i,i are independent of 7. Note that this also implies that the 
leading ©(/"^Z") scaling corrections get suppressed with increasing 7 for any a, including 
a -)■ 1. 

Let us now study the spatial dependence of the entanglement entropy in the 
presence of the trap. For this purpose we define the function 

(19) 



[X 



C-i lim [S^iL/2 -x-L)~ S^{L/2; L)] 



% To get an idea of the accuracy of the universality check, we mention that, using the six available data 
in the range 10^ < I < 10'* to fix the six parameters of the polynomial Ra + ^j=i Cj-e^^^/", we obtain 
i?2 = 0.98892, 0.98892, 0.98892, 0.98889 and R3 = 0.9867, 0.9867, 0.9868, 0.9869 for 7 = 4, 2, 1, 1/2 
respectively. 
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Figure 3. si{x), cf. (fT9l) . vs. X^, for 7 = 1, p = 2, / = 4000 and a few values of the 
size of the chain L. 



This limit is finite and depends only on x and the trap size /. We expect that in the 
TSS limit 



\x 



fa{X), 



X = x/i 



1; 



(20) 

i.e., they are scaling functions of X. Note that the distance of the site L/2 — x from the 
center of the trap is actually y = x + 1/2, but the effects of this difference disappears in 
the TSS limit ^1 —t- 00 at fixed X = y/^i, as 0(^]"^). The large-p limits of fa coincide, 
indeed, since ^1 — )■ L, fITT]) implies 



lim fa{X) = lncos(7rX) 



(21) 



with \X\ < 1/2. Corrections to scaling are again expected to be 0(/~^/"). 

In figure [3] we show data for p = 2, 7 = 1, / = 4000 and a few values of the size of 
the chain L. They show the effects of the finite size L of the chain, and that accurate 
results for larger and larger X require larger and larger L. However, such effects can be 
carefully controlled, so that the results that we will present can be effectively considered 
as infinite-chain results. □ We used chains of size up to L = 4000 for 7 = 1 and up to 
L = 6000 for 7 = 4, which ensure negligible finite-L effects at the trap size considered. 

The numerical results for the von Neumann entropy, reported in figure H] for p = 2 
and p = 4, show that Si(x) approaches a scaling function of X = x/^i in the TSS hmit 
(more precisely, a function of X"^, since it is an even function of X), which is universal, 
i.e., independent of 7. In figure Hlwe also show extrapolations to / — )■ 00 assuming 0(/~^) 
corrections. They are consistent with a universal TSS limit independent of 7. With 

^ The comparison of data in Fig. [3] shows that we can obtain very accurate results around the middle 
of the trap even without requiring I -^ L, a,s one might naively expect. This is essentially related to 
the fact that the relevant length scale of the critical correlations behaves as ^ ~ Z^ and 6* < 1, so that 
we can have ^ ^ L even when I fa L. Note that, with increasing the power p of the confining potential, 
9 increases and approaches one, thus larger chain sizes are required at fixed trap size (we recall that 
for p — >■ 00 we have the formal correspondence L — 21). 
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increasing p, the curves approach the p — )■ 00 hmit 
fast at small X, where we find 



211) . The approach is particularly 



h{X) = c,X' + 0{X') (22) 

with C2 ~ —4.988 for p = 2 and C2 ~ —4.939 for p = 4 (see figure [5]), to be compared 
with C2 = — 7r^/2 ^ —4.93480 for p — )■ cxd, obtained by expanding equation fl2Tl) . 

Analogous results are obtained from the Renyi entropies, see figures |6] and [TJ 
respectively for a = 2 and a = 3. In these cases the approach to the large-Z limit 
is slower, because corrections are 0(/~^/"). As a consequence the I — )■ 00 extrapolations 
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universality of the TSS. For comparison we also show the p — > cx) limit ()2ip . 
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Figure 7. The function 53(2;), cf. ([T9|) . vs. X'^, for p = 2, 7 = 1, 4, and a few values 
of /. Curves for / — >• 00 are obtained by using the data for the largest values of I, 
extrapolating at fixed X using the polynomials X]n=o Cnl~^^^^ with n = 2 for 7 = 4 
and h — 3 for 7 = 1. For comparison we also show the p ~> 00 limit pip . 



tend to be more imprecise with increasing a, see indeed the case a = 3 in figure [71 Larger 
values of a would require larger trap sizes to get reliable / — )■ 00 extrapolations, due to 
the slowly decaying 0(/~^/") scaling corrections. Comparing the results for a = 1, 2, 3, 
we see that the functions /q(X) turn out to depend little on a. 



3. The ID hard-core Bose-Hubbard model with a trapping potential 

We now consider the ID Bose-Hubbard (BH) model ([I]), which is of experimental 
relevance because it describes quasi-one- dimensional cold atoms in optical lattices, see, 
e.g., Refs. [HI ini [ini im [E]- We consider its hard-core limit f/ — )■ 00, which implies 
that the particle number is restricted to the values rij = 0, 1, and allows us to study the 
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effects of the confining potential by very accurate numerical results. 

The hard-core limit of the ID BH model can be mapped into the XX chain model 
with a space-dependent transverse external field, 

^xx = -^ E (^^'^^+1 + ^f^i) - ^ Et-" + y^^^M^ (23) 

i i 

(which is the Hamiltonian ([5]) with 7 = 0, apart from a factor 2). The Pauli spin 
matrices are related to the boson operators of the BH Hamiltonian ifeH; cf. ([T]), by 
'^f = H + hi, erf = i{h\ — hi), erf = 1 — 2h\hi. One can then map the XX chain into 
a model of free spinless fermions by a Jordan- Wigner transformation, given by (jS]) for 
7 = (apart from a factor 2). 

In the absence of the trap, the ID hard-core BH model has three phases: two Mott 
insulator phases, for /i > 1 with {ui) = (empty state) and for /^ < — 1 with (ui) = 1, 
separated by a gapless superfiuid phase for |/i| < 1 characterized by the filling factor 

/ = {ui) = (l/7r)arccos/i. (24) 

A general analysis of the trap-size dependence and TSS at the Mott transitions and 
within the superfiuid phase of the hard-core BH model ([1]) has been presented in 
Ref. [22|. Here we focus on the behaviour of the quantum entanglement within the 
superfiuid phase, i.e., for \fi\ < 1, whose continuum limit is described by a free massless 
bosonic field theory with dynamic exponent -2 = 1, thus a conformal field theory with 
0=1. Specifically, we consider the model at /x = 0, corresponding to half filling in the 
absence of the trap. 

3.1. Modulated multi-TSS in the superfiuid phase 

The trap-size dependence shows subtle effects in the parameter region where the 
homogeneous model without trap has a nonzero filling /, i.e., for yU < 1, including 
the superfiuid phase. This is essentially related to the presence of level crossings of the 
lowest states at finite trap size. They arise because the particle number is conserved, i.e., 
the particle number operator A^ = Yli f^i commutes with the BH Hamiltonian ([T]) even 
in the presence of the trapping potential; thus the eigenvectors do not depend on //, even 
though the eigenvalues do. In the presence of the trapping potential ([2]), the particle 
number A^ = (A^) is finite and increases as A^ ~ / with increasing the trap size /, see fl29|) . 
Therefore, as / — )■ 00, there is an infinite number of ground-state level crossings where 
A^ jumps by 1 and the gap vanishes. In spite of the presence of these level crossings, the 
length scale of the critical modes diverges only the in the large trap-size limit. Note that 
the hard-core limit, f/ — )■ 00 in ([T]), or the spatial dimension do not play any special role, 
so we expect that level crossings at finite trap size are a general feature of the BH model 
([T]) in the presence of a confining potential, when the homogeneous limit of infinite trap 
size has a finite particle density. The main effect of the infinite level crossings in the 
limit / — )■ 00 is that the asymptotic power law behaviours get modulated by periodic 
functions of the trap size /, giving rise to a modulated TSS. 
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Moreover, in the superfluid region, the scahng of the spatial dependence of the 
correlation functions is characterized by two independent length scales having different 
power-law divergences in the large trap-size limit. One of them scales as 

e. ~ I (25) 

and describes the behaviour of observables related to smooth modes, such as the gap 
and the half-lattice entanglement entropy; the other one scales as 

^F-^S C = P/(P + 1), (26) 

and it is found in observables involving the modes at the Fermi scale kp = vr/, such 
as the density-density correlation. This gives rise to a multi-TSS behaviour, which is 
more complicated than that described by the scaling Ansatz ([7]) which applies to the 
XY chain model around its quantum critical point. 

As we shall see, the behaviour of the quantum entanglement of parts of the system 
in its ground state reflects the above general features. 

In the following we briefly mention a number of results obtained in Ref. [22] which 
will be useful for the discussion of the quantum entanglement. 

The particle density of the ID hard-core BH model approaches its local density 
approximation (LDA) in the large-Z limit, i.e., the value of the particle density of the 
homogeneous system at the effective chemical potential 

fi,s{x)=fi+{x/ir. (27) 

The LDA of the particle density reads 

{0 for /icff(a;) > 1, 

(l/7r)arccoSyUeff(a;) for -1 < fics{x) < 1, (28) 
1 for /icff(a;) < -1. 

Corrections are generally suppressed by powers of the trap size, and present a nontrivial 
scaling behaviour [22] in term of the scaling variable Y = x/l^. 

Asymptotically, the total particle number is obtained by integrating the local 
density approximation of the particle density pida: 



N={N) = c^l + 0{l), 


c^ = 2 pida(y) dy. 


(29) 


At /i = we have c^=o ~ 0.76276 for p ■- 


= 2, c^=o ~ 0.859407 for p = 


= 4, and c^=o — )■ 1 


for p -> oo. 







In the superfluid region the interval between subsequent level crossings /q (where 
k enumerates them with increasing trap size) tends to a constant value in the large 
trap-size limit, given by [22] 

hm (/^'^ - 4'0 = Vc,, (30) 

with 0(/~^) corrections (this formula has been checked with great accuracy by numerical 
computations). In the limit p — )■ cxo the interval of level crossings converges toward the 
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Figure 8. The asymptotic beliaviour of tlie gap and the particle density at the 
center of the trap at fi — 0, thus / — 1/2. In formulae [12], A — 7Tt{ip)/^e and 
(^o) — / = (1 — 'P)/^e, where ^e is the entanglement length scale defined in ((38|) . These 
results apply for any p, including p — )■ oo. 



corresponding value in the finite-size behaviour of the homogeneous system with open 
boundary conditions, which is 1/(2/). 

The amphtudes of the asymptotic power law behaviours turn out to be modulated 
by periodic functions of the trap size, whose period is related to the level crossing interval 
f pOj) . In order to describe such phenomenon, since some quantities distinguish between 
odd and even level crossings, we introduce the phase-like variable 



/-/, 



V 



{2k) 




,(2fc+2) _ A2k) 



for /r^</</ 



(2fc+2) 
; 



(31) 



thus < v? < 2. if provides the normalized distance from the largest even level crossing 
smaller than the given /o Values in the range < </) < 1 corresponds to ground states 
with odd particle numbers, while the range 1 < v^ < 2 corresponds to even particle 
numbers. Note that ( !30ll implies that y? is a periodic function of / in the asymptotic 
regime. 

In the superfluid phase, i.e. |/i| < 1, the asymptotic behaviours of the gap and the 
particle density at the center of the trap are given by 

A ~ t{^)l-\ t(^) = 1/2 - 111 - ^1 - 1/2|, 

K)-/~(i-^)r\ 

respectively. The amplitudes of the asymptotic power-law behaviours get modulated by 
periodic functions of the trap size, through the dependence of ip of their amplitudes. 
The values ip = 1/2 and ip = 3/2 correspond to the odd and even peaks of the gap, 
where the particle number A^ = (N) is odd and even respectively. In figure [S] we show 
their behaviours at /i = 0. 

An analogous modulation of the asymptotic behaviours is found in the limit p — )■ oo, 
i.e. for a homogeneous chain of size L with open boundary conditions [22]. In this case. 



In Ref. |22] the variable ip was denoted by 
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Tlie dashed lines show the linear behaviour ±37r/(2^e) for V' = 1 T 1/2, in agreement 
with (I37l). 



the variable y? reads 

^^2{[(L + l)/ + l]/2}, (33) 

where {x} = x — [xj is the fractional part of x (i.e., the sawtooth function). The 
asymptotic behaviour of the gap is given by 



A = 7r(l 



t^'f'X^)L^'- 



(34) 



At /i = 0, i.e., / = 1/2, the variable if takes only discrete values: v^ = 0, 1/2, 1, 3/2 (the 
values 1/2 and 3/2 corresponds to even chain size, with odd and even particle number 
respectively). 



3.2. Asymptotic trap-size dependence of the bipartite entanglement entropies 

Let us consider a chain with even L sites and open boundary conditions, and a trap 
of size / centered between the middle sites of the chain. We again divide it into two 
subsystems and study the behaviours of the von Neumann and Renyi entropies of one of 
the subsystems. We present results up to trap sizes / = O(IO^), corresponding to total 
particle numbers of O(IO^), see equation ( l29l) . We used chains of size up to L = 3000 
to achieve negligible finite-L effects (a more detailed discussion of this point for the XX 
chain can be found in the App. A of Ref. [22]). 

In the absence of the trap, the behaviour of the entanglement entropies is given by 
fITTD with c = 1, with 



Inv^ 



(35) 



^2 + ln(4/7r) + |/„, 

obtained using the results of Refs. |H |28l [29]. Corrections are again expected to be 
0(L~^/°) [3ll[30]. The large-L behaviour of numerical results is indeed compatible with 
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Figure 10. The function Ao{tp) as obtained from data at / w 998, normalized so that 
Ao{l/2) ^ 1 and Ao(3/2) = -1. 



the following asymptotic behaviour 



SaiL/2-x;L) = C. 



ln(L + 1) + lncos(7rX) + e^ 



+ 0(L- 



2/aN 



(36) 



{LcosTiXy/^ 

with X = x/{L + 1), and h = 37r/2 ^ 4.71239, 62 ~ 4.79256 and 63 ^ 4.19726. Note 
that the factor L + 1 in the logarithm and in the definition of X is essential to get 
a formula which is correct up to 0(L~^) at fixed X in the case of the von Neumann 
entropy (a = 1). 

The half-lattice von Neumann entropy in the presence of the trap was already 
considered and numerically investigated in Ref. j22], finding an entanglement length 
scale which increases linearly with the trap size in the whole superfiuid region. 
Specifically, at ;U = and for any power p of the potential, including p — )• 00, the 
half-lattice von Neumann entropy behaves as 



•^i, 



1/2 



lnee + ei + A(¥,)— + 0(C') 



where 



ee 



Og'; 



^^l> 



pT 



l + p 
2p 



(37) 



(38) 



M^) 



^e provides an entanglement length scale, and Ao{(f) is a periodic function of / which 
can be parametrized by 

l + g(y.-l/2) for < y. < 1, 

-1 - g(3/2 -if) for 1 < ^ < 2 ^ ' 

Note that Ao{l =F 1/2) = ±1. In figure M "we show results for p = 2, 4, 00 and 
(p = 1/2, 3/2. They clearly support equation fl37j) . The function Ao{(p) is shown in 
figure [in] for p = 2 for trap sizes around / ~ 998. We estimate q ~ —0.152 for p = 2 
and q ~ —0.044 for p = 4. The parametrization ( 139|) of Ao{<f) is convenient when 
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-1/q 



for (^ = 1 T 1/2, 



comparing with the p — t- 00 hmit, because (/9 = 1/2, 3/2 are the only possible values for 
the homogeneous system in a chain with even L, corresponding to odd and even L/2 
(thus odd and even particle number A^) respectively. The behaviour of the homogeneous 
system with open boundary condition, cf. ( TTTl) . is obtained by replacing ^e — ^ -^ + 1 at 
the values ip = 1/2, 3/2. 

We now show that the half-lattice Renyi entropies show analogous behaviours. In 
figure [11] we show numerical results of the quantity 

H^ = ej''{C-'Sa,,/2 - In^e - e«) (40) 

for p = 2, 4, 00 and ip = 1/2, 3/2. The data appears to follow curves which are 
apparently independent of p, but clearly depending on ip. Therefore, analogously to 
the case of the von Neumann entropy, the asymptotic behaviour of the Renyi entropies 
show the same form for any value of the potential power p, including p — t- 00, when they 
are expressed in terms of the entanglement length scale ^e, i-e., 

Sa,i/2 = a [In^e + e, + &aA,,o(<^)C'/" + 0(C'/")] (41) 

where ba are the same constants entering the behaviour of the half-lattice Renyi 
entanglement entropies of the homogeneous system, cf. (J36l) . and Aa^o is given by (!39l) . 
with very small values of q, i.e. q ~ —0.003 for a = 2, p = 2 and q ~ —0.001 for a = 3, 
p = 2. Thus Aa,o ~ (-1)^ for a > 1. Equations ([37]) and (jH]) imply that the length 
scales defined in IHM behave as 

^a = ^e [1 + 0(C'/°)] (42) 

for any a. 

We now analyze the spatial dependence of the bipartite entanglement entropies. 
We again consider the functions defined in flT9|l . In figures [121 [131 and [Hj we show 
results for p = 2, 4, a = 1, 2, 3, ip = 1/2, 3/2, and several values of the trap size. In 
all cases the behaviours of the functions Sq,(x), defined as in (TT9l) . are characterized by a 
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Figure 12. Sa vs. X^, for a = 1, 2, 3, p = 2, ip — 1/2, and a few values of I. The 
vertical dotted lines show the value X = Xc = l/a^, cf. 



smooth function /q,(X) with X = x/^g, and oscillations around it, which appear larger 
and larger with increasing a, but get suppressed with increasing ^e- These numerical 
results suggest 

s.(x) = /„(X) + 0(C'/"), X^x/^e, (43) 

where the power- law of the corrections can be guessed from the fact that 0(L~^/") 
corrections are expected in homogeneous system with open boundary conditions [31l |30] . 
Moreover, the numerical results strongly suggest that a unique function describes the 
leading behaviour for any a, i.e., /q(X) = f{X) is independent of a, for any p. 
This becomes more evident when considering the average Sa,a of the functions Sa{x) 
at subsequent odd and even peaks of the gap (indeed, as we shall see below, the curves 
at (f = 1/2 and ip = 3/2 have opposite oscillating corrections at least for sufficiently 
small values of X). With increasing p the function f{X) appears to tend to the p — )■ oo 
limit given by (1211) . Actually, the small-X behaviour at p = 2 is already very close to 
the p — 7- oo limit (the large-/ extrapolation of the coefficient C2 of the O(X^) term is 
compatible with its p — )■ oo limit C2 = — 7r^/2 already at p = 2, 4). On the other hand, 
oscillations are strongly dependent on a. 

An interesting numerical relation concerning the leading behaviour f{X) is 
illustrated in figure [ID for p = 2, 4, 6, we find that the difference between f{X) 
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and its p — 7- oo limit given by CFT, i.e. linip^oo /(-^) = IncosTtX, is a function of 
Aieff = i\x\/lf = {aeX)P only, i.e., 

/(X)^lncos7rX + A(/ieff) (44) 

where the (p-independent) function A behaves as A{fics) = 0(/ipff) at small X (more 



with a 



precisely the numerical results indicate A{^es) = a/Xg^ + 0(/igg) with a ~ —1.0). 

Note that the von Neumann and Renyi entropies rapidly decrease around X = Xc 
corresponding to Xc = I, thus Xc = l/oe, cf. fl38|l : Xc is the value of x where 
fJ^esix) = {x/iy = 1, i.e., the value of the chemical potential corresponding to the 
superfluid to empty state transition, where the particle density of the ground state 
vanishes, see equation ( l28ll . We thus expect that, for generic values of /i and p, the 
region around x = Xc develops critical modes related to a low-density Mott transition 
from the superfluid phase to the empty state. The effective chemical potential can be 
expanded around Xc as 



^ieS=ix/iy=l+p- 



l 



0[{x 






(45) 



Thus, the behaviour around Xc is essentially analogous to that arising at the low-density 
Mott transition /x = 1 in the presence of a linear potential Vi ~ r/l. Around Xc, critical 
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{x/l)P for 



modes should appear with length scale ^ ^ l'^ , where a is the exponent associated with 
a linear external potential. The value of a can be inferred by RG arguments analogous 
to those leading to the determination of the trap exponent 6 at the low-density Mott 
transition [20], which give a = 1/3|^ We thus expect that the transition region around 
X = Xc enlarges as 

Ax ~ /^/^ (46) 

independently of the power-law p of the confining potential. This behaviour is confirmed 
by the numerical values of Ax = / — x^ax? where x^ax is the abscissa of the rightmost 
maximum of Sa{x), interpolated to continuous x; figure HE] shows the perfect agreement 

tt The exponent a can be determined by a RG analysis of the perturbation corresponding to a linear 
potential Vi{x) = ux, i.e., Jd'^xdtVi{x)\(l){x,t)\'^, at the fixed point of the continuous theory describing 
the Mott transition |19) . The exponent a is related to the RG dimension 2/„ of the parameter u, which 
can be obtained from the relations y„ — 1 = d + z — y\^\2 = y^ = 2, thus j/„ = 3, and therefore 
a = l/yu = 1/3 for d = 1 and d = 2. 
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with die]). 

The behaviours of the subleading and oscillating 0(^^^/") corrections are instead 
more subtle. In order to study them, we consider the function 

Qa{x)=ej''[s^{x)-si{x)] (47) 

for a > 1. Assuming that the leading behaviour is given by the same function f{X) for 
any a, Qa should provide the 0(^e ) term in ( l43l) . The functions Qaix) for a = 2, 3, 
p = 2, 4, and (p = 1/2, 3/2 are shown in figures IT7 \ [T8] and [T9 | as obtained by numerical 
results at fixed trap size for 10^ < / < 10'^. They clearly show the expected even-odd 
spatial effect proportional to (— 1)"^' = e^''^^'^, but also a (/j-dependent scaling limit as a 
function oi Y = x/l'^ with Q = p/{p + 1). More precisely, the numerical results show 
that the asymptotic behaviour of the the functions Qa{x) can be written as 

QM - {-iTgaAY^ V) + 9a,s{Y, ^) (48) 

apart from further suppressed scaling corrections. Note that this strong evidence of 
scaling provides also an indirect and accurate check of the initial hypothesis that the 



Scaling of entanglement in trapped one- dimensional lattice systems 



21 



-10, 




-5- 



Y Y 

Figure 18. The function Q2{x), cf. (gT]), vs. Y = xl-^l^ , for p = 4 and ip = 1/2, 3/2. 



e, 



-10 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 /«998 - 


- / = 455 - 


■ /«248| 




, 


'11 


even x ~^ 

_ 9 = 1/2 

odd X 







1 


10 


-1 /«999 -- 


■ /«456 - 


- /«249| 






\ / 


AAArtft 




odd X ^""^^ 


5 


- 


Y 


Yn 




. (p=3/2 


A 


A A 





even x ^^_^ 


1 


1 < < < < 




{ 


) 


1 


2 3 



Y Y 
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leading behaviour as a function of X does not depend on a, as in the p — ;■ oo hmit. 
Summarizing, we have 

s.{x) = f{x) + c'/" [{-iTgaAy, V) + g^AY. v)] + 0(c'/") (49) 

where X = x/C,e and Y = x/l^ . We note that the behaviour of the von Neumann 
entropy, a — )■ 1, is quite similar to that of the particle density, which turns out to 
be [22] 

K) ^ pida(X) + c' {{-iTKiY, v) + hsiY, ^)] (50) 

Moreover, as shown by figures [20] and [211 for p = 2 and p = 4 respectively, the 
oscillations appear to be very close to periodic functions when plotted versus V^^^ . 
This is particularly verified in the case of the harmonic potential, where the data at 
the odd and even peaks of the gap, corresponding to ground states with odd and even 
particle number A^, accurately fit the simple function 

Qa{x)\^=l/2,,/2 ~ (-l)'^&a [(-l)^'cOs(cy3) _ l] (51) 

where h^ are the same constants entering the behaviour of the Renyi entropies of the 
homogeneous system, and c ~ 2/3. In figure [221 we show data for other values of (p. 
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Figure 20. The functions Qa{x), cf. (gT]), vs. V^ = x^/P, for p = 2, (p = 1/2, and 
a = 2, 3. 




Figure 21. The functions Qa(a;), cf. gT]), vs. F^ ee xV^'^, for p = 4, (p 
a = 2, 3. 



1/2, and 



Apart from the prefactor (—1)^, there is a small residual dependence on ^j, which is 
apparently relegated to a small dependence of c on 79. 

4. Conclusions 



In this paper we have investigated the scaling behaviour of the bipartite entanglement 
of ID lattice systems in the presence of a power-law confining potential V{r) = {r/iy, 
at a quantum critical point. We have considered lattice systems whose quantum critical 
behaviours are described by 2D conformal field theories in the absence of the trap, thus 
showing logarithmically divergent entanglement entropies, see, e.g., (TTT]) . 

We have considered ID lattice models of even size L and open boundary conditions, 
with a trap of size / centered between the middle sites of the chain. We have divided the 
chain in two parts of length Ia = L/2 — x and L — I a, and computed their von Neumann 
and Renyi entropies Sa{L/2 — x; L). In the presence of a trap of size /, the dependence 
on L disappears when L — )■ 00, leaving only the dependence on x and the trap size /. 
We have studied the scaling behaviour of the bipartite entanglement with increasing the 
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Figure 22. The functions Qa{x), cf. (gll), vs. F^ = x^/l'^, for p = 2, a = 2, 3, several 
values of Lp, for sites with even and odd parity, and trap sizes I sa 998. Qa as a function 
of if has jumps for integer i^; we denote linii^x^„ hy tp = ri^ and lun^y^n hy (p = n~ . 



trap size, keeping the other model parameters fixed at the quantum critical point of the 
homogeneous system without trap. 

As a theoretical laboratory, we have considered the quantum XY chain in an 
external space-dependent transverse field, acting as a trap for the spinless fermions 
of its quadratic Hamiltonian representation. This model presents a quantum critical 
transition belonging to the 2D Ising universality class, thus described by a CFT with 
central charge c = 1/2. In the presence of the trap and for large trap sizes, the quantum 
critical behaviour can be described in the framework of the TSS theory [20], where the 
quantum critical behaviour is cast in the form a TSS, cf. ([7]), with a nontrivial trap 
exponent 9 = p/{p+l), which determines how the length scale of the critical modes at 
the critical point diverges with increasing trap size, i.e., ^ ~ l^ . Exploiting the quadratic 
spinless fermion representation of the XY chain, we have computed the bipartite von 
Neumann and Renyi entanglement entropies up to / ~ 10^. Our results show that at 
the quantum critical point they behave as 

"in^i + ln/2„ + e, + /„(X) + 0(er'^" 



^, = a^^"n' [I + o(/-' 



a 



x = x/ei, wei = ^a + o(er^'") 



-1/qn 



(52) 
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where Ca = c(l + a~^)/12, c = 1/2 is the central charge, and e^ is the same constant 
entering the entanglement entropies of homogeneous systems with open boundary 
conditions, cf. ( IT3l) . The entanglement length scales ^a are defined by imposing 
Sa{L/2;L) = Cq (In^Q + Cq,); the length scales ^a, and in particular ^i derived from 
the von Neumann entropy, behave consistently with the predictions of the TSS theory, 
with 9 = p/{p+l) and Oi independent of 7. The asymptotic ratios Ra and the functions 
/a(X) depend on the power of the potential, show a small dependence on a, but are 
universal with respect to the parameter 7, see figures HJ E], and [71 Moreover, in the 
p ^ 00 limit they tend to fa{X) = IncosvrX with X = x/L, which is the behaviour in 
the homogeneous system with open boundary conditions. 

We have studied the scaling of the entanglement of confined particle systems 
described by the ID BH model ([T]) in the presence of a power-law confining potential, 
with the center of the trap in the superfiuid phase. This model is of experimental 
relevance because it describes cold atoms in quasi-one- dimensional optical lattices. We 
have computed bipartite entanglement in the hard-core limit of the model, which allows 
very accurate numerical calculations by exploiting a map into a lattice model of free 
spinless fermions. In the absence of the trap, the continuum theory corresponding to 
the gapless superfiuid phase is a 2D massless boson field theory, thus a CFT with c = 1. 
In this region the trap-size dependence presents new features with respect to the XY 
chain: a multiscaling behaviour with the existence of different length scales which diverge 
with distinct power laws in the TSS; the presence of level crossings of the lowest states 
at finite trap size, which gives rise to periodic modulations of the amplitudes of the 
asymptotic power-law behaviours, with a period given by the asymptotically constant 
interval between the level crossings. We have shown that the bipartite entanglement 
entropy presents the same features as well. We have computed the von Neumann and 
Renyi entropies at /i = corresponding to half filling / = 1/2 in the absence of the 
trap, up to trap sizes / ~ 10'^, corresponding to total particle numbers A^ = 0(10'^). Our 
results provide a robust numerical evidence of the following scaling behaviour 

\imL^^Sa{L/2 -x;L) = C^yn^e + e^ + f{X) 

+ C'/" [(-i)^'w(r,^) + ^„,,(y,^)] + o(r'/")}, (53) 

where e^ and Oe are constant given in equations (!35l) and (138|) respectively. The leading 
spatial dependence f{X) turns out to be independent of a; it depends on p, approaching 
the CFT dependence f{X) = lncos(7rX) with X = x/L when increasing p, which is 
the result holding in a homogeneous chain of size L. The 0(^e ) term presents an 
even-odd spatial effect proportional to (—1)^ = e^^'^^^, scaling as a function of F = x/l'' 
with ( = p/{p+ 1), modulated by the variable ip defined in equation (13T|) . see figures \T7\ 
[T8l [T9I For the harmonic potential, an accurate description of 0(^e ) term is given 
by (-1)^6« [(-l)^'cos(cy3) _ 1], see figure [201 

The above scaling behaviours have been mostly obtained by analyzing the scaling of 
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numerical (practically exact) results at finite trap size. Of course, it would be worthwhile 
to derive analytically the asymptotic trap-size dependence of the entanglement entropy, 
which may lead to further physical insights. 
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Abstract. 

We study the effects of a power-law trapping potential on the scaling behaviour of 
the entanglement at the quantum critical point of one-dimensional (ID) lattice particle 
systems. We compute bipartite von Neumann and Renyi entropies in the presence of 
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a limited spatial region of the optical lattice created by laser-induced standing waves. 
The theoretical framework [18] is based on the Bose-Hubbard (BH) model [19] in the 
presence of a confining potential coupled to the particle density, i.e., 

H^n = --^(6l6,+6]fe.)+ -5^ni(ni-l)+/i^ni + J]\/(r,)ni, (1) 

(ij) i i i 

where (ij) is the set of nearest-neighbour sites, rii = h\hi is the particle density operator. 
We consider a power-law trapping potential 

V{r) = yPrP = {r/lf, (2) 

where r = \x\, v and p are positive constants and / = 1/v is the trap size. Experiments 
are usually set up with a harmonic potential, i.e., p = 2. Far from the origin the potential 
V{r) diverges, therefore (uj) vanishes and the particles are trapped. The inhomogeneity 
due to the trapping potential strongly affects the phenomenology of quantum transitions 
in homogeneous systems. The correlation functions develop critical behaviours only in 
the limit of large trap size. This critical behaviour can be described in the framework 
of the trap-size scaling (TSS) theory [20, 21]. 

In this paper we focus on the quantum entanglement in lattice particle systems 
confined by an external potential, such as the BH model (1). We investigate how the 
scaling behaviour of the entanglement entropies determined by CFT changes in the 
presence of the inhomogeneity induced by the confining potential. 

As a theoretical laboratory we consider the quantum XY chain in an external space- 
dependent transverse field, acting as a trap for the spinless fermions of its quadratic 
Hamiltonian representation. This model presents a quantum critical behaviour 
belonging to the 2D Ising universality class, thus described by a CFT with central 
charge c = 1/2. In the presence of the trap and for large trap sizes, the quantum critical 
behaviour shows TSS [20, 21]. The presence of the trap induces a length scale C^ at the 
critical point, which scales nontrivially with increasing the trap size, i.e., as .^ ~ /^ with 
a trap exponent 6 = p/{p + 1). 

We then consider the ID BH model (1) in its hard-core limit, i.e., U -^ oo, and study 
the bipartite entanglement entropy in the gapless superfiuid phase, whose continuum 
limit is described by a 2D massless bosonic field theory, thus a CFT with c = 1. In 
this case the trap-size dependence is made more complicated by new features [22]: (i) 
a multiscaling phenomenon, i.e., the existence of different length scales which diverge 
with distinct power laws in the TSS; (ii) the presence of level crossings of the lowest 
states at finite trap size, which gives rise to a peculiar modulation of the amplitudes of 
the asymptotic power-law behaviours. 

We analyze the scaling behaviour of bipartite von Neumann and Renyi 
entanglement entropies in these ID lattice models. We consider chains with even sites 
L and open boundary conditions, and a trap of size / centered between the middle sites 
of the chain. We divide the chain in two connected parts of length I a and L — l^, and 
consider their von Neumann entropy 

Si{Ia,L) = Si{L-Ia;L) = -Ti[pA\npA] (3) 
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and their Renyi entropies 

S^{Ia\L) = S^{L-Ia^L) = ^—InTrp^ (4) 

I — a 

where pA is the reduced density matrix of one of the two subsystems. Renyi entropies are 

useful because they provide information on the spectrum of pa [23]. The von Neumann 

entropy (3) is obtained by taking the hmit a — * 1 in (4). In the presence of a trap of 

size /, the bipartite entanglement entropies of subsystems of size I a = L/2 — x have a 

finite L ^ oo limit, depending on x and /. We study their scaling behaviour in the 

large-/ limit. Note that the p -^ oo limit of the confining potential becomes equivalent 

to a homogeneous chain of size L = 21 with open boundary conditions (more precisely, 

to a chain with even L = 2[/J sites). Thus, in the p -^ oo limit we must recover the 

entanglement entropies of homogeneous systems with open boundary conditions, which 

are determined by CFT. 

In our numerical study we exploit the quadratic spinless fermion representations of 
the XY chain and the ID hard-core BH model, which allow us to perform computations 
for very large systems, since they only require the diagonalization of a L x L matrix 
where L is the number of lattice sites. We compute the entanglement entropies using the 
method outlined in Refs. [24, 25], which exploits the quadratic fermionic representation 
and can be also applied in the presence of the trapping potential. We obtain numerical 
results for chains of size L and open boundary conditions, with a trap of size / centered 
between the middle sites of the chain; we choose L large enough to have negligible 
finite-L effects; we are able to obtain very accurate results for / up to O(IO^). 

Our results show that the bipartite entanglement entropies at the critical point 
diverge logarithmically with increasing the trap size. Moreover, they present notable 
scaling behaviours in the trap-size scaling limit. 

We mention that the bipartite entanglement entropies in the XY and XX model with 
gradients, i.e., in the presence of a linear external field, have been recently investigated 
in Ref. [26]. 

The paper is organized as follows. In Sec. 2 we presents results for the quantum 
XY chain in an external space-dependent transverse field. Sec. 3 is dedicated to the ID 
hard-core BH model. In Sec. 4 we draw our conclusions. 

2. Bipartite entanglement in the XY chain with a space-dependent 
transverse field 

The quantum XY chain in a transverse field is a standard theoretical laboratory for 
issues related to quantum transitions, see, e.g., Ref. [27]. An inhomogeneity analogous 
to the one arising from a trapping potential in particle systems can be achieved by 
considering a space-dependent transverse field, i.e.. 



H- 



XY 



E ^i(i + ^)^«'^^>i + (1 - ^)^'<i] - ^^^ - ^(^«)^^ (5) 
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where a* are the Pauh matrices, 7 7^ and V{x) = {\x\/iy, cf. (2). The XY chain 
can be mapped into a quadratic Hamiltonian of spinless fermions by a Jordan- Wigner 
transformation, obtaining 



H=^ ][clA,c, + -{clB,,c] + h.c.)], 



Aij = 26ij - 5j+ij - Sij+i + 2Q{xi)6ij, /g\ 

Bij = —7 (Si+ij — 5jj+i) , 

Q{x) = Ji + V{x), fi = fi — 1. 

In this picture fi plays the role of chemical potential for the fermion c-particles, and the 
space-dependent field V{x) acts as a trap for the c-particles, making their local density 
[ui) = (qCj) vanish at large distance. In the following, / will be considered as the trap 
size. 

In the absence of the trap, the model undergoes a quantum transition at /i = 
/i — 1 = in the 2D Ising universality class, separating a quantum paramagnetic phase 
for /i > from a quantum ferromagnetic phase for /i < 0. Around /I = 0, the quantum 
critical behaviour shows a diverging length scale, ^ ~ \lA~'^ ■> ^"^^ ^ vanishing energy 
scale, A ~ ^~^, where z and u are universal critical exponents: z = 1 and u = l/y^i = 1 
{y^ is the RG dimension of /i) . 

The presence of the external field V{x) gives rise to a space inhomogeneity, thus 
affecting the critical behaviour of the homogeneous system, which is only recovered in 
the limit of large trap size. In this limit the critical behaviour can be described in the 
framework of the trap-size scaling (TSS) theory [20]. The TSS Ansatz for the asymptotic 
behaviour of the free-energy density in the presence of a confining potential (2) is 

F(/i, T, /, x) = /-'(^+^)J^(/i/'/^ Tl'\ xl-'), (7) 

where x is the distance from the middle of the trap, and 6 is the trap-size exponent which 
determines how the length scale of the critical modes at the critical point diverges with 
increasing trap size, i.e., ^ ~ l^. RG arguments and analytic calculations show that [20] 

e = p/{p + i). (8) 

The TSS limit can be analytically derived within the quadratic spinless fermion 
representation (6), obtaining a Schrodinger-like equation for the lowest states, after 
introducing the rescaled spatial variable X = 7~^/p/~^x, and the rescaled deviation 
from the critical point fir = 7~^/^/i (see Ref. [20] for details). The 7 dependence gets 
absorbed by such rescalings leading to a universal TSS limit. Therefore at the quantum 
critical point /i = and T = 0, any length scale related to the critical modes must have 
the asymptotic behaviour 

^ = aY^^l' [1 + 0(/-^)] (9) 

where the amplitude a depends on the power p of the potential, but not on 7. The 
power law of the scaling corrections is inferred from the analysis of the corrections to 
the continuum Schrodinger-like equation describing the TSS limit. 
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We want to study the quantum entanglement of the ground-state in the presence 
of the external space-dependent potential (2). For this purpose we analyze the trap-size 
dependence of the bipartite von Neumann and Renyi entanglement entropies, defined 
in equations (3) and (4), in the framework of the TSS theory. 

Before presenting the results in the presence of the trap, let us recall some 
established results for the homogeneous XY chain with open boundary conditions, which 
represents the formal p — > oo limit of the model with the confining potential. They can 
be obtained by exploiting the CFT that describes its continuum limit. Setting 

Ia = L/2-x (10) 

(we consider a chain with an even number of sites L) , the Renyi entropies at the critical 
point /i = behave as [4] 



So,{L/2-x;L) ^C, 
where 



In L + In cos ( ^- ) + e^ 



:iii 



"^dt 



C^« = ci±|^ (12) 

and c is the central charge: c = 1/2 for the 2D Ising universality class. The constant Cq, 
can be derived using the results of Refs. [4, 28, 29, 20], obtaining 

Ca = ln7 + ln(8/7r) + y^, 

' 6 / 1 _ _±_\ _J_ _ ^_ J ^ (13) 

1 — q;~2 \asmht/a sinhty sinht 
The behaviour of the von Neumann entropy can be obtained by taking the limit a ^ 1 
of (11). Corrections to the asymptotic behaviour (11) are expected to be 0(L~^/") 
at fixed x/L [30, 31]. They become logarithmic for a -^ oo. In the case of the von 
Neumann entropy, i.e., a — *> 1, corrections are 0(L^^).| 

In the presence of the trap of size /, the dependence on L rapidly disappears for 
L sufficiently larger than the trap size /. We define entanglement length scales C,a from 
the half-lattice entanglement entropies. We write (at /i = 0) 

5^,1/2 = lim Sa{L/2; L) = C^ (In^. + e«) , (14) 

where Cq is the constant defined in (13). We identify the argument of the logarithm as 
the entanglement length scale ^q,, which will generally depend on a and on the power p 
of the potential (2). However, by construction, their p ^ oo limit is ^q, — ;► 2/ = L for 
any a. Since the entanglement length scales are clearly related to the critical modes, 
they are expected to have the asymptotic behaviour (9) at the critical point, although 
their scaling corrections may turn out to be different, depending on a. 

The approach to their asymptotic behaviour C,a ~ l^ niay be guessed by recalling 
that the p ^ oo limit is equivalent to a homogeneous chain of size L = 21 with open 
boundary conditions, where the corrections to the asymptotic behaviour (11) of the 

I It is worth mentioning that, using the relation [29] Sa(lA',L)\y=i ~ ^Sa{'2lA','2L)\~^=o and results 
from Refs. [32, 22], we obtain Si{L/2;L) = ^ [inL + ei + (2 - 37r)/(4L) + 0(1/L2)] for 7 = 1. 
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Figure 1. The ratio S.i/i'j^^^l^) vs ^"^ for p = 2, 4, several values of 7, and trap sizes 
in the range 10^ < / < lO"*. The lines are drawn to guide the eyes. The crosses along 
the y-axis represent the 7-independent I — > 00 extrapolations, see text. 



Renyi entanglement entropies are [31] 0(L~^/"). Then, the correspondence L ^ ^ r^ l^ 
suggests the presence of 0(^~^/"), thus 0(/~^/°'), corrections. In the case of the length 
scale derived from the von Neumann entropy, i.e., for a — i> 1, this argument leads to the 
same scaling corrections of equation (9), which are expected for any length scale defined 
from standard correlation functions, such as Gij = {a^a^). 

The scaling (9) is fully supported by numerical results for the half-lattice von 
Neumann entropy. § Figure 1 shows the results of ^1 for p = 2 and p = 4, several 
values of 7, and for trap sizes in the range 10^ < / < 10^. In order to estimate the 
amplitude ai of the asymptotic behaviour 

6 = ai7'/"/^[l + feiir^ + ...], (15) 

we fit the data of the ratio ^i/{'j^^^l^) to the Ansatz ai + 6i_i/^^ + 61,2^^^^? which turns 
out to be optimal, and check the stability of the results by increasing the minimum 
value of / allowed in the fit. The analysis of the p = 2 data gives ai ~ 1.97431 
for all 7 = 6,4,2,1,1/2,11 confirming its independence of 7, and 611 = &i,i7^^^ with 
bii ^ —2.677 independent of 7. For p = 4 we obtain ai ^ 2.18316 for both 7 = 4, 1. 

In the case of the entanglement length scales C,a defined from the Renyi entropies, 
we generally expect ^q, ~ l^ asymptotically, with 0(/~^/") corrections. Assuming 
universality of the TSS limit, we expect that the ratios ^a/^i approach universal values 

U^l = Ra + 0(/-^/") (16) 

§ The effects of the finite size L of the chain are under complete control in our numerical calculations 
at fixed trap size. This is easily achieved by comparing data at fixed trap size for increasing chain 
size, due to the fast convergence to the infinite-chain limit, see also later. We choose L large enough 
to have negligible finite-L effects, so that the results that we present can be effectively considered as 
infinite-chain results with great accuracy. Therefore, the systematic error related to the finite size of 
the chain is totally negligible in our analyses. 

II Here and in the following the uncertainty on the numerical estimates is at most one on the last 
reported figure. 
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Figure 2. Ratios $a/^i of entanglement lengtli scales vs. I ^/"a in the presence of an 
external potential with p — 2, A, for a = 2, 3, several values of 7, and trap sizes in the 
range 10^ < Z < 10''. 



and therefore 



ic 



y/^/^i + 6M(7r^/" + &«,2(7)r^'/" + •••] , 



(17) 



with aa dependent on p but not on 7. Note that Ra ^ 1 and a^ — i> 2 for p ^ 00. 

This scenario is confirmed by numerical results for the ratios ^q/Ci, see figure 2 
where we show results for p = 2, 4 and a = 2, 3. They fit the general behaviour 

Uei = i?a + $^c,ri^/", (18) 

J=l 

and provide consistent values of Ra for different values of 7, supporting universality. 
For example, we find R2 ^ 0.9889, R^ ^ 0.987 for p = 2, and R2 ^ 0.9970 and 
i?3 ~ 0.996 for p = 4, which are slightly smaller than one, and appear to approach one 
with increasing p.*^ 

Moreover, the results for p = 2 indicate that the amplitudes ha^i of the leading 
0(/^^/") scaling corrections in equation (17) behave as ha,i = &a,i7~''^" where n = 4/3 
and ba,i is independent of 7, with great numerical accuracy. Analogous results are found 
for p = 4, but with k = 6/5. Taking into account that ba,i ~ 7"^/°" for p — ;> 00 (as 
inferred by numerical computations), we may guess n = {p + 2)/{p + 1). Therefore, we 
have that the ratios b^, i/&i,i are independent of 7. Note that this also implies that the 
leading 0(/~^/") scaling corrections get suppressed with increasing 7 for any a, including 
a —>■ 1. 

Let us now study the spatial dependence of the entanglement entropy in the 
presence of the trap. For this purpose we define the function 



(x) = C-' lim [Sail/ 2 



x-L)-S^{L/2-L)] 



(19) 



^ To get an idea of the accuracy of the universality check, we mention that, using the six available data 
in the range 10^ < / < 10^ to fix the six parameters of the polynomial Ra + X^^i Cje"-'^/", we obtain 
R2 = 0.98892, 0.98892, 0.98892, 0.98889 and R^ = 0.9867, 0.9867, 0.9868, 0.9869 for 7 = 4, 2, 1, 1/2 
respectively. 
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/ = 4000, L = 4000 
/ = 4000, L = 3000 
/ = 4000, L = 2000 




Figure 3. si(x), cf. (19), vs. X^, for 7 = 1, p == 2, ^ = 4000 and a few values of the 
size of the chain L. 



This limit is finite and depends only on x and the trap size /. We expect that in the 
TSS limit 



\x 



UX), X^x/i 



1-, 



(20) 



i.e., they are scaling functions of X. Note that the distance of the site L/2 — x from the 
center of the trap is actually y = x + 1/2, but the effects of this difference disappears in 
the TSS limit ^1 — i> 00 at fixed X = y/^i, as 0(^f ^). The large-p limits of /q, coincide, 
indeed, since ^1 ^ L, (11) implies 



lim fa{X) = lncos(7rX) 



(21) 



with \X\ < 1/2. Corrections to scaling are again expected to be 0(/~^'"). 

In figure 3 we show data for p = 2, 7 = 1, / = 4000 and a few values of the size of 
the chain L. They show the effects of the finite size L of the chain, and that accurate 
results for larger and larger X require larger and larger L. However, such effects can be 
carefully controlled, so that the results that we will present can be effectively considered 
as infinite-chain results. "•" We used chains of size up to L = 4000 for 7 = 1 and up to 
L = 6000 for 7 = 4, which ensure negligible finite-L effects at the trap size considered. 

The numerical results for the von Neumann entropy, reported in figure 4 for p = 2 
and p = 4:, show that si(x) approaches a scaling function of X = x/C,i in the TSS limit 
(more precisely, a function of X"^, since it is an even function of X), which is universal, 
i.e., independent of 7. In figure 4 we also show extrapolations to / — ;> cx) assuming 0(/~^) 
corrections. They are consistent with a universal TSS limit independent of 7. With 

'^ The comparison of data in Fig. 3 shows that we can obtain very accurate results around the middle 
of the trap even without requiring I <^ L, as one might naively expect. This is essentially related to 
the fact that the relevant length scale of the critical correlations behaves as £, ^ l^ and 6* < 1, so that 
we can have £, -^ L even when I ra L. Note that, with increasing the power p of the confining potential, 
9 increases and approaches one, thus larger chain sizes are required at fixed trap size (we recall that 
for p — » 00 we have the formal correspondence L = 21). 
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Figure 4. The function si(a;), cf. (19), vs. X"^ , for p = 2, 4, 7 = 1, 4, and a few values 
of I. Curves for / -^ 00 are obtained by using the data for the three largest values 
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universality of the TSS. For comparison we also show the p —f 00 limit (21). 
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The overall agreement with a universal value of C2 is very accurate, leading to the 
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increasing p, the curves approach the p 
fast at small X, where we find 



h{X)=C2X^ + 0{X' 



cxD limit (21). The approach is particularly 



(22) 



with C2 
with C2 



-4.988 for p = 2 and C2 
-7rV2 ^ -4.93480 for p - 



—4.939 for p = 4 (see figure 5), to be compared 

00, obtained by expanding equation (21). 

Analogous results are obtained from the Renyi entropies, see figures 6 and 7, 

respectively for a = 2 and a = 3. In these cases the approach to the large-/ limit 

is slower, because corrections are 0(/~^/°). As a consequence the / — ^ 00 extrapolations 
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tend to be more imprecise with increasing a, see indeed the case a = 3 in figure 7. Larger 
values of a would require larger trap sizes to get reliable I — * 00 extrapolations, due to 
the slowly decaying 0(/^^/") scaling corrections. Comparing the results for a = 1, 2, 3, 
we see that the functions /q,(X) turn out to depend little on a. 



3. The ID hard-core Bose-Hubbard model with a trapping potential 

We now consider the ID Bose-Hubbard (BH) model (1), which is of experimental 
relevance because it describes quasi-one- dimensional cold atoms in optical lattices, see, 
e.g., Refs. [6, 9, 10, 11, 17]. We consider its hard-core limit U — * 00, which implies 
that the particle number is restricted to the values nj = 0, 1, and allows us to study the 
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effects of the confining potential by very accurate numerical results. 

The hard-core limit of the ID BH model can be mapped into the XX chain model 
with a space-dependent transverse external field, 



H- 



J 



XX 



E (^^^^^+1 + ^'<i) - ^ El-" + ^(^^)]^^ (23) 



(which is the Hamiltonian (5) with 7 = 0, apart from a factor 2). The Pauli spin 
matrices are related to the boson operators of the BH Hamiltonian -f^BH, cf. (1), by 
'^i = ^i ~^ bi, <^i = ^ibl — hi), af = 1 — 26|6j. One can then map the XX chain into 
a model of free spinless fermions by a Jordan- Wigner transformation, given by (6) for 
7 = (apart from a factor 2). 

In the absence of the trap, the ID hard-core BH model has three phases: two Mott 
insulator phases, for /x > 1 with (nj) = (empty state) and for fj, < —1 with {rii) = 1, 
separated by a gap less superfluid phase for |/i| < 1 characterized by the filhng factor 

/ = (rii) = (l/7r)arccoS;U. (24) 

A general analysis of the trap-size dependence and TSS at the Mott transitions and 
within the superfluid phase of the hard-core BH model (1) has been presented in 
Ref. [22]. Here we focus on the behaviour of the quantum entanglement within the 
superfluid phase, i.e., for |/i| < 1, whose continuum limit is described by a free massless 
bosonic field theory with dynamic exponent z = 1, thus a conformal field theory with 
c = 1. Specifically, we consider the model at /i = 0, corresponding to half filling in the 
absence of the trap. 

3.1. Modulated multi-TSS in the superfluid phase 

The trap-size dependence shows subtle effects in the parameter region where the 
homogeneous model without trap has a nonzero filling /, i.e., for /i < 1, including 
the superfluid phase. This is essentially related to the presence of level crossings of the 
lowest states at finite trap size. They arise because the particle number is conserved, i.e., 
the particle number operator N = Yli'^i commutes with the BH Hamiltonian (1) even 
in the presence of the trapping potential; thus the eigenvectors do not depend on /x, even 
though the eigenvalues do. In the presence of the trapping potential (2), the particle 
number A^ = (A^) is finite and increases as A^ ~ / with increasing the trap size /, see (29). 
Therefore, as / ^ 00, there is an infinite number of ground-state level crossings where 
A^ jumps by 1 and the gap vanishes. In spite of the presence of these level crossings, the 
length scale of the critical modes diverges only the in the large trap-size limit. Note that 
the hard-core limit, f/ ^ 00 in (1), or the spatial dimension do not play any special role, 
so we expect that level crossings at finite trap size are a general feature of the BH model 
(1) in the presence of a confining potential, when the homogeneous limit of infinite trap 
size has a finite particle density. The main effect of the infinite level crossings in the 
limit / ^ 00 is that the asymptotic power law behaviours get modulated by periodic 
functions of the trap size /, giving rise to a modulated TSS. 
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Moreover, in the superfluid region, the scahng of the spatial dependence of the 
correlation functions is characterized by two independent length scales having different 
power-law divergences in the large trap-size limit. One of them scales as 

is ~ / (25) 

and describes the behaviour of observables related to smooth modes, such as the gap 
and the half-lattice entanglement entropy; the other one scales as 

e^^/S C=P/(P+1), (26) 

and it is found in observables involving the modes at the Fermi scale kp = vr/, such 
as the density-density correlation. This gives rise to a multi-TSS behaviour, which is 
more complicated than that described by the scaling Ansatz (7) which applies to the 
XY chain model around its quantum critical point. 

As we shall see, the behaviour of the quantum entanglement of parts of the system 
in its ground state reflects the above general features. 

In the following we briefly mention a number of results obtained in Ref. [22] which 
will be useful for the discussion of the quantum entanglement. 

The particle density of the ID hard-core BH model approaches its local density 
approximation (LDA) in the large-/ limit, i.e., the value of the particle density of the 
homogeneous system at the effective chemical potential 

fXe{i{x)=fi+{x/lY. (27) 

The LDA of the particle density reads 

{0 for /icff(a;) > 1, 

(l/7r)arccos;Ucff(a;) for -1 < fi^six) < 1, (28) 
1 for ficsix) < -1. 

Corrections are generally suppressed by powers of the trap size, and present a nontrivial 
scaling behaviour [22] in term of the scaling variable Y = x/l^ . 

Asymptotically, the total particle number is obtained by integrating the local 
density approximation of the particle density pida^ 

N^{N) = c^l + 0(1), c^ = 2 / pida(y) dy. (29) 

Jo 

At // = we have c^=o ~ 0.76276 for p = 2, c^=o ~ 0.859407 for p = 4, and c^=o ^^ 1 

for p -^ oo. 

In the superfluid region the interval between subsequent level crossings Iq (where 

k enumerates them with increasing trap size) tends to a constant value in the large 

trap-size limit, given by [22] 

hm f/^'^ - 4'0 = Vc., (30) 



fc— >oo 

with 0(/~^) corrections (this formula has been checked with great accuracy by numerical 
computations). In the limit p ^ oo the interval of level crossings converges toward the 
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corresponding value in the finite-size behaviour of the homogeneous system with open 
boundary conditions, which is 1/(2/). 

The amphtudes of the asymptotic power law behaviours turn out to be modulated 
by periodic functions of the trap size, whose period is related to the level crossing interval 
(30). In order to describe such phenomenon, since some quantities distinguish between 
odd and even level crossings, we introduce the phase-like variable 



/-/ 



^ 



{2k) 




/, 



(2fc+2) j{2k) 



for /r^</</ 



(2fc+2) 



/ 



(31) 



^ <- ^ f-O 

thus < if < 2. if provides the normalized distance from the largest even level crossing 
smaller than the given l* Values in the range < y? < 1 corresponds to ground states 
with odd particle numbers, while the range 1 < v^ < 2 corresponds to even particle 
numbers. Note that (30) implies that </) is a periodic function of / in the asymptotic 
regime. 

In the superfluid phase, i.e. |/i| < 1, the asymptotic behaviours of the gap and the 
particle density at the center of the trap are given by 

(32) 



t{^)l-\ t(<^) = l/2-||l-v9|-l/2|, 

(no) - / ~ (1 - v^)r\ 

respectively. The amplitudes of the asymptotic power-law behaviours get modulated by 
periodic functions of the trap size, through the dependence of y? of their amplitudes. 
The values ip = 1/2 and (p = 3/2 correspond to the odd and even peaks of the gap, 
where the particle number A^ = (A^) is odd and even respectively. In figure 8 we show 
their behaviours at /i = 0. 

An analogous modulation of the asymptotic behaviours is found in the limit p ^ oo, 
i.e. for a homogeneous chain of size L with open boundary conditions [22]. In this case, 

* In Ref. [22] the variable (p was denoted by 4>. 
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the variable y? reads 

^^2{[(L + l)/+l]/2}, (33) 

where {x} = x — [xj is the fractional part of x (i.e., the sawtooth function). The 
asymptotic behaviour of the gap is given by 



A = 7r(l 



/i 



V/^t(^)L- 



(34) 



At /.i = 0, i.e., / = 1/2, the variable </? takes only discrete values: y? = 0, 1/2, 1, 3/2 (the 
values 1/2 and 3/2 corresponds to even chain size, with odd and even particle number 
respectively) . 



3.2. Asymptotic trap-size dependence of the bipartite entanglement entropies 

Let us consider a chain with even L sites and open boundary conditions, and a trap 
of size / centered between the middle sites of the chain. We again divide it into two 
subsystems and study the behaviours of the von Neumann and Renyi entropies of one of 
the subsystems. We present results up to trap sizes / = O(IO^), corresponding to total 
particle numbers of O(IO^), see equation (29). We used chains of size up to L = 3000 
to achieve negligible finite-L effects (a more detailed discussion of this point for the XX 
chain can be found in the App. A of Ref. [22]). 

In the absence of the trap, the behaviour of the entanglement entropies is given by 
(11) with c = 1, with 



e„ = In a/1 -yu2 + ln(4/7r) + y„. 



(35) 



obtained using the results of Refs. [4, 28, 29]. Corrections are again expected to be 
0(L^^/") [31, 30]. The large-L behaviour of numerical results is indeed compatible with 
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the following asymptotic behaviour 



S^iL/2-x;L)=C, 



-(-UL/2-x_ 



ln(L + 1) + lncos(7rX) + e„ 



0(L 



-2/a^ 



(36) 



{LcosTiXy/^ 

with X = x/(L + 1), and h = 37i/2 ^ 4.71239, 62 ^ 4.79256 and 63 ^ 4.19726. Note 
that the factor L + 1 in the logarithm and in the definition of X is essential to get 
a formula which is correct up to 0(L~^) at fixed X in the case of the von Neumann 
entropy (a = 1). 

The half-lattice von Neumann entropy in the presence of the trap was already 
considered and numerically investigated in Ref. [22], finding an entanglement length 
scale which increases linearly with the trap size in the whole superfiuid region. 
Specifically, at ;U = and for any power p of the potential, including p ^ 00, the 
half-lattice von Neumann entropy behaves as [22] 



5'l,l/2 



lnee + ei + ^(¥.)— + 0(C') 



where 



^e 



a^i, 



v^rl- 



pT 



1+p 
2p 



(37) 



(3^ 



M^) 



^e provides an entanglement length scale, and Ao{if) is a periodic function of / which 
can be parametrized by 

l + g(y.-l/2) for < y. < 1, 

-1 - g(3/2 ~ip) for 1 < y? < 2 ^ ^ 

Note that Ao{l =F 1/2) = ±1. In figure 9 we show results for p = 2, 4, 00 and 
if = 1/2, 3/2. They clearly support equation (37). The function Ao{(p) is shown in 
figure 10 for p = 2 for trap sizes around / ~ 998. We estimate q ~ —0.152 for p = 2 
and q ~ —0.044 for p = 4. The parametrization (39) of Ao{<f) is convenient when 
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p = 2, 4, oo. Tlie dashed lines show the linear behaviour ±6q,^( 
with 62 = 4.79256 and 63 = 4.19726, cf. (41). 



-l/a 



for (^ = 1 T 1/2, 



comparing with the p ^ 00 hinit, because (p = 1/2, 3/2 are the only possible values for 
the homogeneous system in a chain with even L, corresponding to odd and even L/2 
(thus odd and even particle number N) respectively. The behaviour of the homogeneous 
system with open boundary condition, cf. (11), is obtained by replacing ^e — ;> L + 1 at 
the values ip = 1/2, 3/2. 

We now show that the half-lattice Renyi entropies show analogous behaviours. In 
figure 11 we show numerical results of the quantity 

Ha = ee^°(C-^S,,i/2 - In^e - e«) (40) 

for p = 2, 4, 00 and ip = 1/2, 3/2. The data appears to follow curves which are 
apparently independent of p, but clearly depending on ip. Therefore, analogously to 
the case of the von Neumann entropy, the asymptotic behaviour of the Renyi entropies 
show the same form for any value of the potential power p, including p — »• 00, when they 
are expressed in terms of the entanglement length scale ^e? i-e., 

Sa,i/2 = C^ [In^e + e, + 6aA,,„(¥p)C'/" + 0(C'/")] (41) 

where ba are the same constants entering the behaviour of the half-lattice Renyi 
entanglement entropies of the homogeneous system, cf. (36), and A^^o is given by (39), 
with very small values of q, i.e. q ~ —0.003 for a = 2, p = 2 and q ~ —0.001 for a = 3, 
p = 2. Thus Aa^o ~ (~1)^ for a > 1- Equations (37) and (41) imply that the length 
scales defined in (14) behave as 

e„ = ^e [1 + 0(C'/")] (42) 

for any a. 

We now analyze the spatial dependence of the bipartite entanglement entropies. 
We again consider the functions defined in (19). In figures 12, 13, and 14 we show 
results for p = 2, 4, a = 1, 2, 3, (p = 1/2, 3/2, and several values of the trap size. In 
all cases the behaviours of the functions Sq,(x), defined as in (19), are characterized by a 
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Figure 12. Sa vs. X^ , for a = 1, 2, 3, p = 2, (p = 1/2, and a few values of I. The 
vertical dotted lines show the value X = Xc = l/oe, cf. (38). 



smooth function fa{X) with X = x/^e, and oscillations around it, which appear larger 
and larger with increasing a, but get suppressed with increasing ^e- These numerical 
results suggest 

s„(x) = /,(X) + 0(C'/"), X^x/^e, (43) 

where the power-law of the corrections can be guessed from the fact that 0(L^^/") 
corrections are expected in homogeneous system with open boundary conditions [31, 30]. 
Moreover, the numerical results strongly suggest that a unique function describes the 
leading behaviour for any a, i.e., fa{,X) = f{X) is independent of a, for any p. 
This becomes more evident when considering the average Sa,a of the functions Sq,(x) 
at subsequent odd and even peaks of the gap (indeed, as we shall see below, the curves 
aX ip = 1/2 and if = 3/2 have opposite oscillating corrections at least for sufficiently 
small values of X). With increasing p the function f{X) appears to tend to the p ^ oo 
limit given by (21). Actually, the small-X behaviour at p = 2 is already very close to 
the p ^ oo limit (the large-/ extrapolation of the coefficient C2 of the O(X^) term is 
compatible with its p ^ oo limit C2 = —tt'^/2 already at p = 2, 4). On the other hand, 
oscillations are strongly dependent on a. 

An interesting numerical relation concerning the leading behaviour f{X) is 
illustrated in figure 15: for p = 2, 4, 6, we find that the difference between f{X) 
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Figure 13. s^ vs. X'^ for a = 1,2,3, p = 2, and / w 10^, aX ^p = 1/2, ^p = 3/2, 
corresponding to A^ = 761, 762, and their average Sa,a = (sq|i^=i/2 + Sa|^=i/3)/2, 
which has much smaller oscillations. 



and its p — !► CX3 limit given by CFT, i.e. linip^oo /(-^) = IncosvrX, is a function of 
/ieff = {\x\/lf = {a.Xy only, i.e., 

/(X) ^ In cos vrX + A(/ieff) (44) 

where the (p-independent) function A behaves as A{jjics) = 0(/^gg) at small X (more 
precisely the numerical results indicate A{fics) = OyU^g + 0(/i^g) with a ~ —1.0). 

Note that the von Neumann and Renyi entropies rapidly decrease around X = Xc 
corresponding to Xc = I, thus Xc = 1/ae, cf. (38); Xc is the value of x where 
fJ'esix) = {x/iy = 1, i.e., the value of the chemical potential corresponding to the 
superfluid to empty state transition, where the particle density of the ground state 
vanishes, see equation (28). We thus expect that, for generic values of /i and p, the 
region around x = Xc develops critical modes related to a low-density Mott transition 
from the superfluid phase to the empty state. The effective chemical potential can be 
expanded around Xc as 



/ieff = (x/iy = l+p- 



l 



0[{x — Xc 



(45) 



Thus, the behaviour around Xc is essentially analogous to that arising at the low-density 
Mott transition /x = 1 in the presence of a linear potential Vi ~ r/l. Around Xc, critical 
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{x/l)P for 



modes should appear with length scale ^ ~ Z'^, where a is the exponent associated with 
a linear external potential. The value of a can be inferred by RG arguments analogous 
to those leading to the determination of the trap exponent 6 at the low-density Mott 
transition [20], which give a = 1/3. jj We thus expect that the transition region around 
X = Xc enlarges as 

Ax ~ /^/^ (46) 

independently of the power-law p of the confining potential. This behaviour is confirmed 
by the numerical values of Ax = I — Xmax, where Xmax is the abscissa of the rightmost 
maximum of Sa{x), interpolated to continuous x; figure 16 shows the perfect agreement 

tt The exponent a can be determined by a RG analysis of the perturbation corresponding to a linear 
potential Vi{x) — ux, i.e., /d'^xdt Vi(a;)|(/)(x, i)p, at the fixed point of the continuous theory describing 
the Mott transition [19]. The exponent a is related to the RG dimension y^ of the parameter u, which 
can be obtained from the relations j/u — 1 = d + z — y\M2 = y^ = 2, thus j/„ = 3, and therefore 
a = l/yu = 1/3 for d = 1 and d = 2. 
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fits with the form al^/^, cf. (46). 
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with (46). 

The behaviours of the subleading and oscillating 0(^^^/") corrections are instead 
more subtle. In order to study them, we consider the function 

Qa{x)=ej''K{x)-s,{x)] (47) 

for a > 1. Assuming that the leading behaviour is given by the same function f{X) for 
any a, Q^ should provide the 0{^e ) term in (43). The functions Qa.{x) for a = 2, 3, 
p = 2, 4, and (p = 1/2, 3/2 are shown in figures 17, 18 and 19, as obtained by numerical 
results at fixed trap size for 10^ < / < 10'^. They clearly show the expected even-odd 
spatial effect proportional to (—1)^ = e^''^^^, but also a c/9-dependent scaling limit as a 
function oi Y = x/l'^ with C = p/(p + !)• More precisely, the numerical results show 
that the asymptotic behaviour of the the functions Qa{x) can be written as 



Qa{x) ^ {-!)'' ga,o{Y^ V^) + 9a,s{Y^V) 



(4^ 



apart from further suppressed scaling corrections. Note that this strong evidence of 
scaling provides also an indirect and accurate check of the initial hypothesis that the 
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leading behaviour as a function of X does not depend on a, as in the p ^ oo hmit. 
Summarizing, we have 

s.{x) = f{x) + c'/" [{-iTgaAy, v) + 9aAy, v)] + 0{C'h (49) 

where X = x/C,e and Y = x/l^ . We note that the behaviour of the von Neumann 
entropy, a ^ 1, is quite similar to that of the particle density, which turns out to 
be [22] 

K) ^ pida(X) + C' [{-^TKiY, if) + KiY, if)] (50) 

Moreover, as shown by figures 20 and 21, for p = 2 and p = 4 respectively, the 
oscillations appear to be very close to periodic functions when plotted versus V^^^ . 
This is particularly verified in the case of the harmonic potential, where the data at 
the odd and even peaks of the gap, corresponding to ground states with odd and even 
particle number A^, accurately fit the simple function 

Qa{x%=l/2,,/2 ~ (-l)^&a [(-1)" COs(cr=^) - l] (51) 

where ha are the same constants entering the behaviour of the Renyi entropies of the 
homogeneous system, and c ~ 2/3. In figure 22 we show data for other values of (/?. 
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Figure 21. The functions Qa{x), cf. (47), vs. Y^ = x^/l'^, for p = 4, y? 
a = 2, 3. 



1/2, and 



Apart from the pref actor (—1)^, there is a small residual dependence on ip, which is 
apparently relegated to a small dependence of c on 79. 

4. Conclusions 



In this paper we have investigated the scaling behaviour of the bipartite entanglement 
of ID lattice systems in the presence of a power-law confining potential V{r) = {r/iy, 
at a quantum critical point. We have considered lattice systems whose quantum critical 
behaviours are described by 2D conformal field theories in the absence of the trap, thus 
showing logarithmically divergent entanglement entropies, see, e.g., (11). 

We have considered ID lattice models of even size L and open boundary conditions, 
with a trap of size / centered between the middle sites of the chain. We have divided the 
chain in two parts of length Ia = L/2 — x and L — Ia, and computed their von Neumann 
and Renyi entropies Sa{L/2 — x; L). In the presence of a trap of size /, the dependence 
on L disappears when L -^ 00, leaving only the dependence on x and the trap size /. 
We have studied the scaling behaviour of the bipartite entanglement with increasing the 
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Figure 22. The functions Qa{x), cf. (47), vs. F^ = a;^/^^, for p = 2, a = 2, 3, several 
values of (f, for sites with even and odd parity, and trap sizes I ~ 998. Qa as a function 
of if has jumps for integer ip; we denote lim;px^„ by (^3 = n+ and lim;^^„ by 99 = n^. 



trap size, keeping the other model parameters fixed at the quantum critical point of the 
homogeneous system without trap. 

As a theoretical laboratory, we have considered the quantum XY chain in an 
external space-dependent transverse field, acting as a trap for the spinless fermions 
of its quadratic Hamiltonian representation. This model presents a quantum critical 
transition belonging to the 2D Ising universality class, thus described by a CFT with 
central charge c = 1/2. In the presence of the trap and for large trap sizes, the quantum 
critical behaviour can be described in the framework of the TSS theory [20], where the 
quantum critical behaviour is cast in the form a TSS, cf. (7), with a nontrivial trap 
exponent 6 = p/{p+ 1), which determines how the length scale of the critical modes at 
the critical point diverges with increasing trap size, i.e., ^ ~ l^. Exploiting the quadratic 
spinless fermion representation of the XY chain, we have computed the bipartite von 
Neumann and Renyi entanglement entropies up to / ~ 10^. Our results show that at 
the quantum critical point they behave as 

In^i + lni?„ + e„ + /„(X) + 0(er'/" 



limL^oo5'„(L/2 - x; L) 

^, = a^^'/n' [I + o(/-' 



Ca 



x = x/i^, U^i = Ra + OiCi 



-1/qx 



(52) 
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where Ca = c(l + a~^)/12, c = 1/2 is the central charge, and e^ is the same constant 
entering the entanglement entropies of homogeneous systems with open boundary 
conditions, cf. (13). The entanglement length scales ^^ are defined by imposing 
Sa{L/2;L) = Cq, (In^Q + Cq); the length scales C,a, and in particular ,^i derived from 
the von Neumann entropy, behave consistently with the predictions of the TSS theory, 
with 6 = p/{p+l) and ai independent of 7. The asymptotic ratios Ra and the functions 
/q,(X) depend on the power of the potential, show a small dependence on a, but are 
universal with respect to the parameter 7, see figures 4, 6, and 7. Moreover, in the 
p ^ 00 limit they tend to fa{X) = IucosttX with X = x/L, which is the behaviour in 
the homogeneous system with open boundary conditions. 

We have studied the scaling of the entanglement of confined particle systems 
described by the ID BH model (1) in the presence of a power-law confining potential, 
with the center of the trap in the superfluid phase. This model is of experimental 
relevance because it describes cold atoms in quasi-one-dimensional optical lattices. We 
have computed bipartite entanglement in the hard-core limit of the model, which allows 
very accurate numerical calculations by exploiting a map into a lattice model of free 
spinless fermions. In the absence of the trap, the continuum theory corresponding to 
the gap less superfluid phase is a 2D massless boson field theory, thus a CFT with c = 1. 
In this region the trap-size dependence presents new features with respect to the XY 
chain: a multiscaling behaviour with the existence of different length scales which diverge 
with distinct power laws in the TSS; the presence of level crossings of the lowest states 
at finite trap size, which gives rise to periodic modulations of the amplitudes of the 
asymptotic power-law behaviours, with a period given by the asymptotically constant 
interval between the level crossings. We have shown that the bipartite entanglement 
entropy presents the same features as well. We have computed the von Neumann and 
Renyi entropies at /i = corresponding to half filling / = 1/2 in the absence of the 
trap, up to trap sizes / ~ 10^, corresponding to total particle numbers A^ = 0(10'^). Our 
results provide a robust numerical evidence of the following scaling behaviour 

\imL^^Sa{L/2 -x-L) = ajln^e + e„ + /(X) 

+c'/" [{-irgc.Ay.v)+9c,s{y^v)] + o{c"")}. (^3) 

ee = ae/, X = x/i,, r = a:/r/(p+i), 

where Cq, and ag are constant given in equations (35) and (38) respectively. The leading 
spatial dependence f{X) turns out to be independent of a; it depends on p, approaching 
the CFT dependence f{X) = lncos(7rX) with X = x/L when increasing p, which is 
the result holding in a homogeneous chain of size L. The 0{^e ) term presents an 
even-odd spatial effect proportional to (—1)^ = e^^^p^ ^ scaling as a function oiY = x/l^ 
with C, = p/{p+ 1), modulated by the variable ip defined in equation (31), see figures 17, 
18, 19. For the harmonic potential, an accurate description of 0{^e ) term is given 
by i-l)^ba [{-ly cos{cY^) - 1], see figure 20. 

The above scaling behaviours have been mostly obtained by analyzing the scaling of 
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numerical (practically exact) results at finite trap size. Of course, it would be worthwhile 
to derive analytically the asymptotic trap-size dependence of the entanglement entropy, 
which may lead to further physical insights. 

Acknowledgments 

Helpful discussions with Pasquale Calabrese are gratefully acknowledged. 

References 

[1] Calabrese P, Cardy J and Doyon B Eds. 2009 Entanglement entropy in extended quantum systems^ 

J. Pfiys. A: Math. Theor. 42 500301 

[2] Holzhey C, Laursen F and Wilczek F 1994 Nucl. Pfiys. B 424 443 

[3] Vidal G, Latorre J I, Rico E and Kitaev A 2003 Ptiys. Rev. Lett. 90 227902 

[4] Calabrese P and Cardy J 2004 J. Stat. Mech. P06002 

[5] Calabrese P and Cardy J 2009 J. Phys. A: Math. Theor. 42 504005 

[6] Bloch I, Dalibard J and Zwerger W 2008 Rev. Mod. Phys. 80 885 

[7] Greiner M, Bloch I, Mandell M O, Hansch T and Esslinger T 2002 Nature 415 39 

[8] Stoferle T, Moritz H, Schorl C, Kohl M and Esslinger T 2004 Phys. Rev. Lett. 92 130403 

[9] Paredes B, Widera A, Murg V, Mandel O, Foiling S, Cirac I, Shlyapnikov G, Hansch R W and 
Bloch I 2004 Nature 429 277 

[10] Kinoshita T, Wenger T and Weiss D S 2004 Science 305, 1125 

[11] Kinoshita T, Wenger T and Weiss D S 2005 Phys. Rev. Lett. 95 190406 

[12] Hadzibabic Z, Kriiger P, Cheneau M, Battelier B and Dalibard J 2006 Nature 441 1118 

[13] Kriiger P, Hadzibabic Z and Dahbard J 2007 Phys. Rev. Lett. 99 040402 

[14] Foiling S, Widera A, Miillcr T, Gerbier F and Bloch I 2006 Phys. Rev. Lett. 97 060403 

[15] Spielman I B, PhiUips W D and Porto J V 2007 Phys. Rev. Lett. 98 080404 

[16] Spielman I B, Phillips W D and Porto J V 2008 Phys. Rev. Lett. 100 120402 

[17] Clement D, Fabbri N, Fallani L, Fort C and Inguscio M 2009 Phys. Rev. Lett. 102 155301 

[18] Jaksch D, Bruder C, Cirac J I , Gardiner C W and ZoUer P 1998 Phys. Rev. Lett. 81 3108 

[19] Fisher M P A, Weichmann P B, Grinstcin G and Fisher D S 1989 Phys. Rev. B 40 546 

[20] Campostrini M and Vicari E 2010 Phys. Rev. A 81 023606 

[21] Campostrini M and Vicari E 2009 Phys. Rev. Lett. 102 240601 

[22] Campostrini M and Vicari E 2010 Phys. Rev. A 81 063614 

[23] Calabrese P and Lefevre A 2008 Phys. Rev. A 78 032329 

[24] Peschel I 2003 J. Phys. A: Math. Gen. 36 L205 

[25] Latorre J I and Riera A 2009 J. Phys. A: Math. Theor. 42 504002 

[26] Eisler V, Igloi F and Peschel I 2009 J. Stat. Mech. P02011 

[27] Sachdcv S 1999 Quantum Phase Transitions (Cambridge Univ. Press) 

[28] Jin B-Q and Korcpin V E 2004 J. Stat. Phys. 116 79 

[29] Igloi F and Juhasz R 2008 Europhys. Lett. 81 57003 

[30] Calabrese P and Cardy J 2010 J. Stat. Mech. P04023 

[31] Calabrese P, Campostrini M, Essler F and Nienhuis B 2010 Phys. Rev. Lett. 104 095701 

[32] Laflorencie N, Sorensen E S, Chang M-S and Affleck I 2006 Phys. Rev. Lett. 96 100603 

[33] Zhou H-Q, Barthcl T, Fiaerestad J and SchoUwock U 2006 Phys. Rev. A 74 050305 



